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ABSTRACT. 

During  the  project,  using  an  approach  based  on  optimal  control  theory,  we  have  first  implemented  a  pulse 
shaper  design  tool.  Using  as  tuning  parameter  the  longitudinal  dependence  of  the  second  order  response  %(2)  of 
the  medium  where  propagation  takes  place,  the  goal  has  been  achieved  by  minimizing  the  distance  between  a 
given  target  pulse  and  the  second  harmonic  at  output  of  a  device  for  second  harmonic  generation.  The  outcome 
of  this  first  theoretical  demonstration  of  pulse  shaping  technique  is  the  function  %(2)(z)  which  is  necessary  to 
obtain  the  desired  shaping  effect.  Starting  from  results  we  have  previously  obtained  [1],  we  have  demonstrated 
here  both  wavelength  tunability  and  arbitrary  pulse  shaping  capability. 

In  this  framework  we  have  also  explored  ultrabroadband  optical  pulse  propagation  in  nonlinear  quadratic  media 
[2,  3]  to  address  engineered  supercontinuum  generation  as  a  valuable  tool  for  spectral  shaping. 

According  to  the  work  plan  described  in  the  research  proposal,  during  the  third  month  of  the  project  activities, 
Matteo  Conforti  and  Costantino  De  Angelis  have  been  visiting  their  research  partners  in  the  United  States, 
namely: 

•  Prof.  Alejandro  B.  Aceves  and  Prof.  Ildar  Gabitov  at  the  Department  of  Mathematics  of  the  Southern 
Methodist  University,  in  Dallas. 

•  Dr.  Michael  Scalora,  AMSRD-WSS  Charles  M.  Bowden  Research  Facility  at  Redstone  Arsenal, 
Huntsville. 

•  Prof.  Triantaphyllos  R.  Akylas,  Department  of  Mechanical  Engineering,  Massachusetts  Institute  of 
Technology,  Cambridge,  Massachusetts. 

During  the  seventh  month  of  the  project  Stefan  Wabnitz  has  been  visiting  Prof.  Alejandro  B.  Aceves  at  the 
Department  of  Mathematics  of  the  Southern  Methodist  University,  in  Dallas. 

The  research  work  carried  on  during  these  visits  has  been  very  fruitful  and  has  been  a  real  boost  for  the  project; 
in  particular  we  have  identified  there  the  need  for  a  better  understanding  of  plasmonic  wave  dynamics,  thereby 
providing  the  theoretical  underpinnings  of  potential  innovative  applications.  To  this  end,  a  combined  analytical 
/numerical  study  of  light  propagation  in  various  structures  composed  of  metals  and  dielectrics  has  been 
pursued. 

In  this  framework  periodic  structures,  consisting  of  metal  and  dielectrics  have  been  investigated  and  we  have 
demonstrated  the  existence  of  new  solitary  wave  solutions  of  a  model  describing  light  propagation  in 
nonuniform  (linearly  and  nonlinearly)  waveguide  arrays.  This  general  model  describes  energy  localization  and 
transport  in  different  physical  settings,  ranging  from  metal-dielectric  (i.e.  plasmonic)  to  photonic  crystal 
waveguides.  The  solitons  exist  for  both  focusing,  defocusing  and  even  for  alternating  focusing-defocusing 
nonlinearity  [4]. 

IN  DETAIL. 

In  the  field  of  modelling  plasmonic  nanostructures  and  engineering  of  their  guiding  properties,  the  work  done 
in  this  project  has  paved  the  way  to  the  management  of  broadband  light  propagation.  Moreover,  all  the 
numerical  tools  already  developed  by  UNIBS  to  analyse  the  modal  properties  of  metal-dielectric  structures  and 
the  field  evolution  in  the  linear  regime  are  to  be  considered  as  one  of  the  result  of  this  research  project.  As  far 
as  diffraction  management  is  concerned,  we  have  considered  propagation  in  a  structure  where  two  different 
plasmonic  devices  are  involved.  The  first  one  is  a  system  composed  of  alternating  metal  (30  nm)  and  dielectric 
(120  nm)  layers. 

In  figure  1  we  show  the  dispersion  relation  and  the  profiles  of  the  two  modes  supported  by  this  directional 
coupler  (to  be  called  DC1).  In  stark  contrast  with  conventional  waveguides,  the  fundamental  mode  of  this 
structure  is  odd.  Moreover  we  can  see  that  the  fundamental  mode  has  one  node  whereas  the  second  one  has  no 
nodes. 


Figure  1 :  a)  Dispersion  relation  of  the  fundamental  (odd)  and  second  order  (even)  mode  of  the  coupler  DC1 .  b) 
Fundamental  (red)  and  second  order  (blue)  mode  @  600nm. 


A  second  well  known  plasmonic  guiding  structure  is  the  metal  nanoparticle  array,  where  the  energy  transport  is 
supplied  by  electromagnetic  resonant  coupling  between  metal  particles  arranged  in  a  linear  chain.  We  have 
studied  double  nanoparticle  chains,  where  the  electromagnetic  energy  is  confined  between  two  linear  chains. 
As  an  example  we  considered  double  chain  waveguides  composed  of  silver  nanoparticles  with  a  radius  of  50 
nm  in  a  dielectric  host  with  a  longitudinal  separation  of  1 10  nm  and  a  distance  between  the  chains  of  150  nm. 
In  figure  2  we  show  the  dispersion  relation  and  the  mode  profiles  of  two  coupled  waveguides  composed  of 
three  nanoparticle  chains  (to  be  called  directional  coupler  DC2).  As  opposed  to  the  previous  directional 
coupler,  the  fundamental  mode  here  is  even  and  the  second  order  mode  is  odd. 


Figure  2:  a)  Dispersion  relation  of  the  coupler  DC2.  Thick  curves,  fundamental  (even)  and  second  order  (odd)  mode. 
Dashed  curves,  higher  order  modes.  Thin  horizontal  line  denotes  a  wavelength  of  600  nm.  b)  Fundamental  and  second 

order  mode  @  600nm. 


We  have  then  considered  the  uniform  array  obtained  using  directional  couplers  DC1  (DC2)  as  basic  building 
blocks.  The  field  evolution  along  the  waveguide  arrays  has  been  simulated  without  any  approximation  by 
solving  Maxwell's  equations  through  a  frequency-domain  finite-element  method,  using  arrays  composed  by  17 
waveguides.  The  central  waveguide  of  the  array  is  excited  with  a  Gaussian  field,  which  spreads  during 
propagation  and  generates  the  typical  diffraction  pattern  observed  also  in  conventional  waveguide  arrays:  two 
outermost  wings  and  a  few  less  intense  peaks  in  the  central  waveguides.  The  same  qualitative  behavior  is 
observed  for  both  arrays  since  the  intensity  evolution  is  not  influenced  by  the  diffraction  sign.  On  the  other 
hand  the  phase  front  curvature  of  the  propagating  field  depends  on  the  diffraction  sign. 

Therefore,  if  we  alternate  arrays  characterized  by  normal  and  anomalous  diffraction,  the  input  field  shape  can 
be  periodically  recovered  as  demonstrated  in  figure  3. 


y  Lum] 

Figure  3:  Time  average  power  flow  along  propagation  in  the  diffraction  managed  device  @  600nm. 

In  the  nonlinear  regime,  we  have  obtained  solitary  wave  solutions  of  a  model  describing  light  propagation  in 
binary  (linearly  and  nonlinearly)  waveguide  arrays.  This  model  describes  energy  localization  and  transport  in 


various  physical  settings,  ranging  from  metal— dielectric  (i.e.  plasmonic)  to  photonic  crystal  waveguides.  The 
solitons  exist  for  focusing,  defocusing  and  even  for  alternating  focusing— defocusing  nonlinearity. 

In  this  project  we  have  considered  a  binary  array  designed  in  such  a  way  that  the  coupling  between  successive 
waveguides  switches  periodically  from  C  to  -C(l+e),  thus  opening  a  gap  centered  at  zero  Bloch  momentum  in 
the  linear  dispersion  relation.  We  consider  also  a  binary  Kerr  nonlinearity  and  we  look  for  self-sustained 
nonlinear  propagation  in  the  form  of  gap  solitons  in  such  a  structure.  Specifically,  extending  previously  derived 
results,  we  obtain  in  the  continuum  limit  exact  analytical  solutions  for  both  stationary  and  "walking"  gap 
solitons  moving  along  the  spatial  coordinate  with  a  tunable  velocity. 

We  have  also  focused  our  attention  on  the  existence  domain  for  bright  gap  solitons.  These  solitons  belong  to  a 
family  with  two  free  parameters:  the  velocity  v  and  the  energy  unbalance. 


Figure  4:  Existence  conditions  for  gap  soliton  solutions  on  transverse  and  longitudinal  wave  number.  Continuous  blue  line 
refers  to  the  dispersion  relation  of  the  linear  problem;  the  red  region  is  the  existence  domain  for  gap  soliton  solutions. 

Such  solitons  display  several  interesting  and  unusual  features,  unique  to  this  type  of  waveguide  structures,  and 
are  possible  even  in  the  case  of  alternating  focusing-defocusing  nonlinearity  [4]. 
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We  obtain  solitary- wave  solutions  of  a  model  describing  light  propagation  in  binary  (linearly  and  nonlinearly) 
waveguide  arrays.  This  model  describes  energy  localization  and  transport  in  various  physical  settings,  ranging 
from  metal-dielectric  (i.e.,  plasmonic)  to  photonic  crystal  waveguides.  The  solitons  exist  for  focusing,  defocusing, 
and  even  for  alternating  focusing-defocusing  nonlinearity. 


DOI:  10. 1103/PhysRevA.83. 043822 

I.  INTRODUCTION 

Discrete  optics  in  coupled  waveguides  has  been  an  area  of 
intense  research  activities  during  the  last  three  decades  (see  [1] 
for  a  recent  review).  Most  efforts  have  been  devoted  to  the 
analysis  of  linear  and  nonlinear  properties  of  uniform  waveg¬ 
uide  arrays,  i.e.,  arrays  composed  of  equally  spaced  identical 
waveguides,  and  both  one-dimensional  and  multidimensional 
configurations  have  been  considered  theoretically  as  well  as 
experimentally  [2-8]. 

On  the  other  hand,  nonuniform  waveguide  arrays  offer  a 
richer  setting  where  engineering  of  the  periodic  structure  can 
provide  further  degrees  of  freedom.  In  this  context,  zigzag 
waveguide  arrays  (i.e.,  the  cascade  of  arrays  characterized  by 
alternating  tilt  angles)  have  been  introduced  to  get  diffraction 
management  [9].  Binary  arrays  composed  of  waveguides 
with  different  wave  numbers  have  been  thoroughly  studied 
(see  [10-13]  and  related  works)  since  they  exhibit  interesting 
features,  such  as  double  refraction,  due  to  their  intrinsic  two- 
band  nature.  Moreover,  binary  arrays  with  different  coupling 
coefficients  have  been  considered,  since  they  might  offer  a 
more  feasible  experimental  framework  in  which  to  exploit  a 
two-band  structure  in  the  linear  and  nonlinear  regimes  [14,15]. 
In  this  instance,  the  use  of  photonic  crystal  waveguides  [16]  or 
waveguides  based  on  plasmonic  confinement  [17,18]  offers  a 
unique  setting  in  which  to  exploit  propagation  in  the  so-called 
alternating  positive  and  negative  coupling  regime  [19,20]. 
Efremidis  et  al.  [20],  in  particular,  studied  nonlinearly  uniform 
arrays  where  the  coupling  coefficients  are  equal  in  modulus 
but  of  opposite  sign.  In  this  case  the  structure  does  not 
possess  a  gap  in  the  linear  spectrum,  and  can  be  reduced  to 
a  uniform  array  by  a  phase  transformation.  Exploiting  this 
transformation,  families  of  discrete  solitons  were  calculated 
starting  from  the  well-known  discrete  Schrodinger  equation. 

In  this  paper  we  consider  a  binary  array  designed  in  such  a 
way  that  the  coupling  between  successive  waveguides  switches 
periodically  from  C  to  —  C(1  +  e),  thus  opening  a  gap  centered 
at  zero  Bloch  momentum  in  the  linear  dispersion  relation. 
We  consider  also  a  binary  Kerr  nonlinearity  and  we  look  for 
self- sustained  nonlinear  propagation  in  the  form  of  gap  solitons 
in  such  a  structure.  Specifically,  extending  previously  derived 
results  [21-30],  we  obtain  in  the  continuum  limit  exact  analyt¬ 
ical  solutions  for  both  stationary  and  “walking”  gap  solitons 
moving  along  the  spatial  coordinate  with  a  tunable  velocity. 
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PACS  number (s):  42.65.Tg,  73.20.Mf,  78.67.Pt 

II.  THEORETICAL  MODEL  AND  SOLITON  SOLUTIONS 

According  to  coupled  mode  theory  and  taking  into  account 
third-order  nonlinearities  in  the  form  of  a  pure  Kerr  effect,  the 
governing  equations  read  as  [10] 

iKz  +  PnK  +  Cn.XE'n_x  +  Cn+lE'n+l  +  Xn\E'n\2E'n  =  0, 

where  E'n  is  the  amplitude  of  the  modal  field  Mn(x)  of  the  nth 
waveguide;  /3n  is  the  propagation  constant  of  each  individual 
waveguide  (/3W  =  /3  +  A/3/2  for  n  even  and  /3n  =  /3  —  A/3/2 
for  n  odd);  Xn ,  the  site-dependent  nonlinear  coefficient,  is 
Y\  (72)  for  n  even  (odd);  and  C„_i,C„+i  are  the  coupling 
coefficients  with  the  ( n  —  l)th  and  the  ( n  +  l)th  waveguides, 
respectively.  In  the  specific  case  of  interest,  Cn_  1  =  C\  and 
Cn+ 1  =  C2  when  n  is  even,  whereas  Cn-\  =  C2  and  Cn+\  = 
C 1  when  n  is  odd.  We  then  perform  the  transformation  E'n  = 
En  exp (i/3z)  and  we  separately  consider  the  mode  amplitudes 
in  the  even  and  odd  waveguides.  Finally,  E2n  =  An  and 
E2n-i  =  Bn  are  governed  by  the  following  two  sets  of  coupled 
equations  with  constant  coefficients: 

AB  9 

iAnz  H — —  Aw  +  C\Bn  +  Bn+ 1  +  7i|Aw|  An  =  0, 

(1) 

A/3  9 

iBnz - —  An-i  +  Ci An  +  Y2\Bn\  Bn  —  0, 

where  C2  has  been  set  equal  to  1 ,  without  loss  of  generality. 

Assuming  Bloch-wave  disturbances,  ( An,Bn )  oc 
exp {i(nkx  +kzz)},  the  linear  dispersion  relation  of  Eqs.  (1) 
reads 

k2  =  ^  +  C2  +  1  +  2Ci  cos kx. 

Note  that  a  band  gap  opens  whenever  A/3  ^  0  and/or  for 
Ci  /  zb  1,  the  band  edges  corresponding  to  the  wave  number 
kx  =  0  for  Ci  <  0  and  kx  —  ix  for  Ci  >  0.  Moreover,  there 
is  numerical  evidence  [10,28]  that  discrete  solitons  can  reside 
inside  this  gap.  We  shall  make  a  comprehensive  analytical 
study  of  stationary  and  moving  gap  solitons  on  the  basis  of  an 
equivalent  continuous  model. 

Specifically,  for  Ci  <  0,  in  the  neighborhood  of  kx  =  0,  we 
use  the  expansions 

An±\ (z)  =  U(x,z)  ±  Ux(x,z)  +  \uxx(x,z)  H - , 

Bn±i(z)  =  w(x,z)  ±  wx(x,z)  +  \wxx(x,z )  H - 
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to  obtain  (as  a  first-order  approximation) 

A  /3  2 

iuz  H - u  +  wx  +  €w  +  yi \u\  u  =  0, 

2  (2) 

A/3  2 

iwz - —w  —  Ux  +  €U  +  y2\w\  w  =  0, 

where  C\  =  —1  +  6.  This  equation  system  also  arises  in  the 
neighborhood  of  kx  =  ix  for  C\  —  1  +  6,  following  a  similar 
expansion  procedure  after  the  change  of  variables  ( An,Bn )  — > 
(-1  )n(An,Bn). 

We  now  look  for  both  stationary  and  walking  self-confined 
solutions  of  the  system  defined  by  Eqs.  (2).  To  this  end,  we 
use  the  following  trial  functions  [23]: 


u(x,z)  =  \[Kxgx{%)  +  iK2g2{%)\cxp(if  cos  Q), 
w(x,z)  =  ^[Kig!($)  -  iK2g2($)]exp(ii/rcos  Q ), 


X  +  vz 

\jf  = 

vx  +  z 

Vi  —  V* 2 

\/ 1  —  V2 

1  +  <A1/4 

k2 

n  -  v 

l-v)  ’ 

~  U  +  V 

with  g  12  two  arbitrary  complex  functions,  —  1  ^  v  ^  1  and 
0^  Q  ^  n .  Although  not  necessary,  for  the  sake  of  simplicity, 
from  now  on  we  set  A/3  =  0  (i.e.,  the  biatomic  nature  of  the 
array  is  left  to  the  coupling  coefficients  only). 

Substitution  of  the  ansatz  (3)  into  Eqs.  (2)  gives  (s  =  y\  + 
Y2,  d  =  yi  -  y2) 

~gi  +  i  cos(Q)gi  +  ieg2+  |2gi  +  2\g2\2gi  -  gjgt) 

K22\g2\2gi  -  2K2\gl\2g2  +  K2g2g*2)  =  0, 
g2  +  i  cos(Q)g2  +  iegi  +  ^(K^gi^gi  +  2\gx\2g2  -  g2g2) 
~(Kt\gi\2gi+2K2\g2\2gl  -  K2g2gl)  =  0. 


These  equations  have  the  invariant  P  =  |gi|2  —  \g2\2',  as  we 
are  interested  in  bright  solitons,  we  set  P  =  0,  so  that  \g\  |2  = 
\g2\2  and  gi,2(£)  =  /(§)exp[i01>2(§)].  Finally,  using  r\  =  f2 
and  /x  =  0\  —  02,  we  get 

.  _  dH  .  _  dH 

^  9/x  ’  ^  dr]  ’ 

H  =  2r](€  cos  /x  +  cos  Q) 

(4) 

j  9  (K\  K*  \ 

-8"  (T  +  f  +2-cos(2/,)j 

~^r]2(K2  +  K2)  sin /*. 

Equations  (4)  represent  a  one-dimensional  (thus  integrable) 
Hamiltonian  system,  and  solitary-wave  solutions  correspond 
to  the  separatrix  trajectories  that  are  homoclinic  to  (i.e., 
emanate  from  and  return  to)  the  unstable  fixed  points  of  (4).  In 
the  following  we  assume  s  >  0,  since  the  results  can  be  easily 
extended  to  the  case  s  <  0  by  the  substitution  /x  ->  /x  +  7T, 
Q  -*  n  -  Q. 

Bright  solitons  emanate  from  the  unstable  fixed  point 
(VOiPo)  =  (0,  ±  arccos[—  cos(2)A])  and  correspond  to  level 
curves  of  the  Hamiltonian  H(rjo,fjio)  =  0.  By  exploiting  H  = 
0,  we  can  derive  the  expression  of  r\  as  a  function  of  /x  from 
the  definition  of  H : 

16(6  cos  /x  +  cos  Q ) 

s\_lt  +  ^2  +  2  —  cos(2/x)j  +  2 d[K\  +  K%)  sin /x 

(5) 

By  inserting  Eq.  (5)  into  /x  =  it  follows  that 


fi  =  — 2(cos  <2  +  6  cos  /x).  (6) 

This  equation  can  be  easily  integrated  to  obtain 


/x(£)  =  —  2arctan 

and,  upon  substitution  in  Eq.  (5), 


g  +  C°S  ^  lanh|  s/f2  -  cos({9)2£ ] 

6  —  COS  Q 


(7) 


m  = 


sech(<$§)2(cos  Q  +  6) 


-4dbiatanh(8^)-\-sb2[l  +  a2  tanh(<5§)2]  —2s 


[1— a2  tanhQSf  )2]2  ’ 
1+a2  tanh(5f  )2 


(8) 


where  a  =  ^+^1+  5  =  ^e2  -cos(g)2,  b\  =  K2  +  K2, 
b2  —  3 — y  +  3. 

Once  we  get  the  solutions  /7(§)  and  /x(§),  we  can  find  the 
fields  u(x,z )  and  w(x,z)  by  substitution  in  (3)  and  solving 

6\  =  cos  2  +  6  cos /x - ^[Ki  +  2  —  cos(2/x)] 

8 

-  |t7  (3^2  +  /T|)  sin/x.  (9) 


The  expression  of  0\  is  rather  cumbersome  except  for  station¬ 
ary  solutions  (v  =  0).  In  fact  by  inserting  (5)  and  (6)  in  (9)  it 
is  straightforward  to  show  that  0\  =  /x / 2. 

We  now  focus  our  attention  on  the  existence  domain 
for  bright  gap  solitons.  These  solitons  belong  to  a  family 
with  two  free  parameters:  the  velocity  v  (—  l^u^l)  and 

2  (arccos  6  <  Q  <  tv  —  arccos  6).  From  Eq.  (5)  we  note, 

however,  that,  as  the  array  parameters  (e,s,d)  are  changed,  the 

amplitude  r]  can  diverge  at  some  points,  entailing  that  some 

( v ,  2)  couples  are  not  allowed. 


043822-2 


ENERGY  LOCALIZATION  AND  TRANSPORT  IN  BINARY  . . . 


PHYSICAL  REVIEW  A  83,  043822  (201 1) 


More  specifically,  for  a  given  Q ,  we  find  that  rj(%)  is 
bounded  for  all  § ,  and  bright  soliton  solutions  can  exist  only 
above  a  critical  velocity  v;  as  long  as  s  >  \d\  this  critical 
velocity  is  0  so  that  all  possible  v  and  Q  values  in  the 
(Q,v)  plane  can  be  attained;  however,  when  \d\  ^  s  a  different 
situation  arises:  for  arccos e  <  Q  <  7r/2,  solutions  only  exist 
above  a  critical  velocity  vcr : 

\vcr |  =  — A2  -d2  +  ~\d\^J2  -  A  (10) 

.S’ 


For  7t/2  ^  Q  ^  7r  —  arccos  6  such  a  critical  velocity  does 
not  exist;  however,  bright  soliton  solutions  are  permitted  only 
for  Q  >  7T  —  arccos(c  VO  with  t  given  by 


t  = 


s2(3  +  /i)  -  (J/2)2  +  y(J/2)4  -  2s2d2  f 2 (\  +  /1) 


2s2 


/1  = 


V  +  ^2 


(ID 


/2 


k2  +  £22. 


This  last  condition  can  be  derived  by  looking  at 
the  phase  plane  (77,/z).  When  Q  =  n  —  arccos^VO, 
another  unstable  fixed  point  exists  characterized  by  [ip  — 
— sgn(J)  arccos[— cos(2)/6)],  r\p  —  8e2yi  —  cos(g)2/62 
sec(2)/[—  \d\(K\  +  K^)  +  2,71  —  cos(g)2/62]  and  having 
Hamiltonian  H(rjp,fjip)  =  0.  In  this  instance  the  separatrix 
trajectory  is  heteroclinic,  connecting  the  points  (^0,/xo) 
and  ( r}p,iJLp ),  and  the  resulting  solution  corresponds  to  a 
kink  soliton.  If  we  increase  Q  above  this  threshold,  the 
trajectories  in  the  phase  plane  become  unbounded,  preventing 
the  existence  of  localized  solutions. 

The  existence  conditions  on  v  and  Q  can  be  easily 
translated  into  conditions  on  the  soliton  transverse  phase 
kx  ~  r>cos(2)/(Vl  —  v2)  and  propagation  constant  kz  = 
cos(<2)/(7T^2),  as  illustrated  in  Fig.  1. 


FIG.  1.  (Color  online)  Existence  conditions  on  kx  and  kz. 
Continuous  line  refers  to  the  dispersion  relation  of  the  discrete 
problem;  dashed  line  shows  the  dispersion  relation  of  the  continuous 
approximation  and  the  filled  region  corresponds  to  the  existence 
domain  of  the  gap  soliton  solutions  (e  =  0.25,  s  =  2,  d  =  2.1). 


(b) 


waveguide  number 


(d) 


FIG.  2.  (a)  Bifurcation  diagram  of  the  Hamiltonian  system 
(v  =  0):  continuous  line  for  stable  centers,  dashed  line  for  unstable 
saddles;  inset  shows  phase  plane  for  Q  =  1.77  (open  circle  indicates 
saddle,  filled  circle  indicates  stable  center);  (b)  field  amplitude  in 
even  (continuous)  and  odd  (dashed)  waveguides  (v  =  0);  (c)  field 
evolution  along  the  array  for  v  =  0;  (d)  field  evolution  along  the 
array  for  v  =  0.5.  In  all  panels,  e  =  0.25,  Q  =  1.77,  5  =  2,  d  =  0. 


III.  EXAMPLES 

In  this  section,  we  discuss  some  specific  examples  of  the 
soliton  solutions  derived  earlier.  We  also  show  the  robustness 
of  our  solutions  in  some  representative  cases  where  we 
consider  propagation  in  different  arrays. 

As  a  first  example,  we  consider  an  array  with  all  the 
waveguides  having  the  same  nonlinear  response  (s  =  2,  d  = 
0).  In  this  case,  bright  soliton  solutions  do  exist  for  arccos  €  < 
Q  ^  7T  —  arccos  €  as  one  can  also  infer  from  the  bifurcation 
diagram  of  Fig.  2(a),  that  shows  the  amplitude  770  of  the 
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FIG.  3.  (a)  Bifurcation  diagram  of  the  Hamiltonian  system 
(v  =  0):  continuous  line  for  stable  centers,  dashed  line  for  unstable 
saddles;  left  (right)  inset  shows  phase  plane  for  Q  =  1.5  (Q  =  1.77) 
(open  circle  indicates  saddle,  filled  circle  indicates  stable  center);  (b) 
field  amplitude  in  even  (continuous)  and  odd(dashed)  waveguides 
(v  =  0);  (c)  field  evolution  along  the  array  for  v  =  0;  (d)  field  evolu¬ 
tion  along  the  array  for  v  =  0.5.  In  all  panels,  e  =  0.25,  Q  =  1.77, 
5  =  2,  d  =  2. 


FIG.  4.  (a)  Bifurcation  diagram  of  the  Hamiltonian  system 
(v  =  0):  continuous  line  for  stable  centers,  dashed  line  for  unstable 
saddles;  left  (right)  inset  shows  phase  plane  for  kink  (bright)  soliton 
(open  circle  indicates  saddle,  filled  circle  indicates  stable  center);  (b) 
field  amplitude  in  even  (continuous)  and  odd  (dashed)  waveguides 
(v  =  0);  (c)  field  evolution  along  the  array  for  v  =  0;  (d)  field 
evolution  along  the  array  for  v  =  0.5.  In  all  panels,  e  =  0.25, 
Q  =  1.77,  5  =  2,  d  =  2.1. 
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fixed  points  of  system  (4)  as  a  function  of  the  parameter 
Q.  In  fact,  in  this  interval  both  the  unstable  saddle  r) o  =  0 
and  a  stable  center  with  r) o  /  0  exist.  From  the  phase  plane 
depicted  in  the  inset  (corresponding  to  Q  =  1.77),  it  is 
evident  that  the  separatix  describing  the  soliton  emanates  from 
(??o,/xo)  =  (0,  arccos[-cos(g)/6]),  turns  around  the  center 
and  returns  to  the  fixed  point  with  vanishing  r].  Note  that 
in  this  situation  (« d  =  0)  the  v  =  0  case  shows  perfect  mirror 
symmetry  between  the  field  in  the  even  and  odd  sites  [see 
Fig.  2(b)];  as  d  increases  this  mirror  symmetry  is  obviously 
lost.  Figures  2(c)  and  2(d)  show  the  propagation  of  a  stationary 
(v  =  0)  and  a  moving  (v  =  0.5)  soliton. 

The  second  example  we  are  considering  in  this  section 
corresponds  to  s  =  2 ,d  =  2  (i.e.,  an  interlaced  linear-nonlinear 
array).  As  one  can  see  from  Fig.  3(a),  for  this  choice  of 
parameters  bright  soliton  solutions  for  v  =  0  exist  only  for 
7r/2  <  Q  <  n  —  arccos6,  because  the  stable  center  does  not 
exist  for  Q  <  n / 2.  In  this  instance,  the  phase  portrait  is 
qualitatively  different  for  Q  greater  or  less  than  7t/2  [right 
and  left  insets  of  Fig.  3(a)].  For  Q  >  n /2  the  phase  portrait 
is  similar  to  the  d  =  0  case,  except  for  the  asymmetry  with 
respect  to  fi.  For  Q  <  7r/2  the  separatrix  emanating  from  the 
saddle  is  not  closed  and  separates  orbits  of  unbounded  motion 
from  periodic  motion.  As  a  consequence  solitons  do  not  exist. 

It  is  remarkable  to  note  that  even  in  the  case  of  interlaced 
focusing-defocusing  nonlinearities  soliton  solutions  still  exist 
as  clearly  demonstrated  in  Fig.  4  for  s  =  2, d  =  2A\  this 
applies  also  to  solutions  walking  along  the  array  as  shown 
in  Fig.  4(d).  Note  also  that,  as  can  be  seen  from  Fig.  4(a),  in 
this  case  we  do  not  have  bright  soliton  solutions  for  Q  <  1 .73; 
however,  as  we  have  already  noted  above,  in  the  presence  of  a 
nonzero  transverse  velocity  we  have  access  to  this  region  of  Q 
values.  This  is  what  we  can  see  in  Fig.  5  where  propagation  in 
an  interlaced  focusing-defocusing  array  is  shown  for  Q  =  1 .72 
and  v  =  0.5;  note  that,  remarkably,  this  last  case  corresponds 
to  a  situation  where  we  do  not  have  bright  soliton  solutions 
with  zero  transverse  velocity. 

Another  interesting  feature  of  the  interlaced  focusing- 
defocusing  case  is  the  existence  of  flat-top  and  kink  solitons, 
due  to  the  presence  of  an  additional  saddle  in  the  bifurcation 
diagram.  It  is  possible  that  the  two  saddles  possess  the  same 


waveguide  number 


FIG.  5.  Field  evolution  along  the  array  for  v  =  0.5;  here  e  = 
0.25,  Q  =  1.72,  s  =  2,d  =  2.1. 
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(a) 


(b) 

FIG.  6.  (a)  Example  of  a  flat-top  (thin  lines)  and  a  kink  (thick 
lines)  soliton.  Here  e  =  0.25,  Q  =  1.786553  604  650  208  for  the  dark 
soliton  ( Q  =  1.786  553  7  for  the  flat-top  soliton),  s  =  2,  d  =  2.5, 
v  =  0.5.  (b)  Field  evolution  along  the  array  for  the  kink  soliton. 


Hamiltonian:  in  this  case  the  heteroclinic  orbit  connecting  the 
two  points  gives  rise  to  a  kink  soliton  [left  inset  of  Fig.  4(a)]. 
As  Q  approaches  the  existence  limit  defined  by  (11),  bright 
solitons  become  wider  and  eventually  take  a  kinklike  shape. 
An  example  of  this  kind  of  solution  is  reported  in  Fig.  6(a). 
Figure  6(b)  shows  the  propagation  of  the  kink  soliton  with 
velocity  v  =  0.5. 


IV.  CONCLUSIONS 

We  have  analyzed  a  model  describing  light  propagation 
in  a  binary  array,  accounting  for  alternating  positive  and 
negative  linear  coupling  as  well  as  nonuniform  nonlinearity. 
This  model  can  be  applied  in  different  physical  settings 
such  as  plasmonic,  Bragg,  and  photonic  crystal  waveguides. 
We  derived  exact  bright  and  kink  soliton  solutions  in  the 
long-wavelength  (i.e.,  continuous)  limit.  Such  solitons  display 
several  interesting  and  unusual  features,  unique  to  this  type 
of  waveguide  structure,  and  are  possible  even  in  the  case  of 
alternating  focusing-defocusing  nonlinearity. 
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1 .  Equipment  purchased . $  7407 

(Server  Gigabyte  Intel  Urbanna:  SUB  SYSTEM  INTEL  SR2600URBRPR  RACK  2U  2  x  750  WATT  -  2  x 
PROCESSORS  INTEL  XEON  SIX  CORE  2.66GHZ  X5650  -  48GB  MEMORY  DDR3  ECC  REG.  (12  x  4GB  )  - 
RISER  CARD  3  x  BUS  F/H  PCI  EXPRESS  2.0  x8  -  CONTROLLER  RAID  ADAPTEC  3405-SGL  -  2  x  HARD 
DISK  WESTERN  DIGITAL  WD1002FBYS  RAID  EDITION  SATA2  7200RPM  1TB  -  GRAPHICS  ON  BOARD 


LCC  PILOT  II  64MB) 

2.  Software  Purchase  (COMSOL  Multiphysics) .  $1733 

3.  Travel .  $  6165 

4.  Publication  Charges .  $1000 

5.  Overheads .  $  6195 


